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Abstract 


The build up of the stress whose relaxation is presumed 
to account for pulsar frequency glitches can be attributed 
to various mechanisms, of which the most efficient involve 
differential rotation of the neutron superfluid in the inner 
layers of the (magnetically braked) solid crust of a rotat¬ 
ing neutron star. In such a case it is usually supposed 
that the stress is attributable to pinning of superfluid vor¬ 
tices to crust nuclei, but it has been suggested that, even 
if the pinning effect is too weak, a comparably large stress 
might still arise just from the deficit of centrifugal buoy¬ 
ancy in the slowed down crust. The present work is a 
re-examination that investigates the way such processes 
may be affected by considerations that were overlooked in 
the previous work notably uncertainties about the “effec¬ 
tive” masses that have to be attributed to the “free” super¬ 
fluid neutrons to allow for their entrainment by the ionic 
crust material. Though restricted to a Newtonian formu¬ 
lation, this analysis distinguishes more carefully than has 
been usual between true velocities, which are contravari- 
antly vectorial, and so called “superfluid velocities” that 
are proportional to momenta, which are essentially cov- 
ectorial, a technicality that is important when more than 
one independent current is involved. The results include a 
Proudman type theorem to the effect that the superfluid 
angular velocity must be constant on slightly deformed 
Taylor cylinders in the force free case, and it is shown how 
to construct a pair of integral constants of the motion that 
determine the solution for the pinned case assuming beta 
equilibrium. 


1 Introduction 

This work is intended as a comparison of the qualitative 
mathematical consequences of hypotheses of various alter¬ 
native kinds (particularly concerning vortex pinning and 
chemical equilibrium) that have been invoked in the con¬ 
struction of differentially rotating neutron star models of 
the type developed to acccount for phenomena observed 
in pulsars. 

It is hoped that this qualitative study will be helpful 
for the preparation of quantitative models designed to im¬ 
prove on those of the simple kind constructed by Prix et 
al f5T)l , for which (as for the analogously constructed mod¬ 
els m in a relativistic framework) separately conserved 
proton and superfluid neutron currents were able to de¬ 
viate from Proudman - Taylor type corotation illustrated 
on figure n without giving rise to the (glitch producing) 
stresses that are actually anticipated. In earlier work m 
such stresses were obtained by including allowance for the 
effect of vortex pinning, but it has since been remarked m 
that deviations from corotation would give rise to stresses 
in any case (with or without pinning) when allowance is 
made for the consideration that the proton and neutron 
numbers will not be separately conserved but will undergo 
transfusion (by beta processes) so as to achieve chemical 
equilibrium. 

The work of Prix m and of Carter et al m used a 
simplified treatment that neglected the mechanism of rel¬ 
ative entrainment between the differentially rotating con¬ 
stituents, an effect which had already been considered in 
earlier work (both in a Newtonian framework m and also 
m in a relativistic treatment) but which had been consid¬ 
ered to be only of secondary relevance as a minor quan¬ 
titative correction. The effect of entrainment (but not 
transfusion) in the liquid core was allowed for in the more 
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recent work of Comer & Andersson m and of Prix et al. 
m, and it has since been recognised m to be particu¬ 
larly important in the crust layers that are the primary 
concern of the present analysis, which will consider the 
(glitch producing) stresses that can be expected to arise 
from the combined effect of entrainment, pinning, trans¬ 
fusion, and also nuclear stratification. 

Explanations of the glitches of the pulsar period P = 
2n /fi observed in rotating neutron stars may be broadly 
classified into two main categories. In what may be de¬ 
scribed as the “deformation” category, the discontinuous 
spin up Sfl is attributed to a sudden change of the geomet¬ 
ric configuration of the matter distribution ESI, whereby 
the relevant moment of inertia is decreased so that - for 
a given angular momentum - the angular velocity must 
increase. While conceivably sufficient to account for cases 
such as the glitches of order Sfl/fl ~ 10 -8 observed in the 
Crab pulsar, such a deformation mechanism is inadequate 
for explaining the frequent occurrence of the larger gliches 
<5f2/f2 ~ 10“ 6 observed in Vela like pulsars, not to mention 
the enormous glitch Sfl/fl ~ 1.6 x 10 -5 recently observed 
in PSR J1806-2125 by Hobbs et al. ESI - 

To account for the frequent occurence of such large 
glitches, it has long been generally recognised HI to be 
necessary to invoke a mechanism belonging to what may 
be referred to as the “transfer” category, involving inter¬ 
action between two dynamically distinct constituents of 
the star. In such a case the discontinuous increase in the 
observed angular velocity is attributed to transfer of 
angular momentum from a more rapidly rotating internal 
constituent to a more slowly rotating constituent that is 
directly coupled to the outer magnetosphere whose pulsa¬ 
tions are directly observed. However it is not so clear what 
actually happens during the discontinous process involved. 
It has been suggested m that the stress due to pinning of 
neutron superfluid vortices to the ionic crust material may 
build up elastic stress to the point where the solid stric¬ 
ture itself breaks down (as in a terrestrial earthquake). In 
another much studied scenario Emm it is postulated that 
before this can happen there will be a crisis of a different 
kind whereby a threshold is reached at which the vortices 
in the crust become collectively unpinned (whereas in the 
fluid core it is expected that vortices in the core will be 
subject to a strong dissipative drag, which will simulates 
pinning but will not be limited by any threshold). Another 
possibility ED E3 EH El is that pinning in the crust may 
be too weak to be effective, but it has been pointed out 
m that differential rotation might in principle by itself 
engender sufficient elastic stress to cause a breakdown of 
the solid structure even in the absence of pinning. 

The main purpose of the present work is not to draw any 
specific conclusions about the relative plausibility of these 
various postulates in particular cases, but to consider, gen¬ 


erally, how such alternatives will affect the predicted out¬ 
come, and how estimates of relevant quantities such as 
subcomponent moments inertia and stresses may need to 
be modified to take account of some effects that tended to 
be neglected in preceding discussions, including stratifica¬ 
tion and particularly the entrainment between the ionic 
crust material and the “free” superfluid neutrons, which 
thereby aquire “effective masses” that can be defined in 
various ways, but that in any case are likely [ITil fPll ITb] to 
be large enough to make an important difference. 

In order to carry out such a revision, and as a secondary 
purpose in its own right, this work employs and develops a 
Newtonian formulation of the covariant kind mm whose 
use would be taken for granted in a relativistic treatment, 
in which care is taken to distinguish between true veloc¬ 
ities, which are contravariantly vectorial, and so called 
“superfluid velocities” that are proportional to momenta, 
which are essentially covectorial. 

Although deeply rooted in the well known principles of 
Lagrangian and Hamiltonian mechanics, the fundamental 
distinction between velocities and their canonical conju¬ 
gates, namely momenta, has commonly been obscured in 
the context of superfluidity as traditionally formulated El 
in a Newtonian framework, for which it has long been cus¬ 
tomary to use the term “velocity” indiscriminately for any¬ 
thing having the physical dimensions thereoff, regardless 
of its actual mathematical role. The distinction is rela¬ 
tively unimportant for a single constituent perfect fluid or 
superfluid such as that of ordinary Helium at zero temper¬ 
ature, or even for a two contituent fluid in which entrain¬ 
ment is absent ESI, but it becomes important when there 
are two or more independently moving fluids with entrain¬ 
ment, which means that the momenta need no longer be 
aligned with their conjugate velocities. 

As well as simplifying the algebra in the work presented 
below, the introduction of the canonically covariant for¬ 
mulation also has the advantage of providing a helpful 
step towards the development of a fully General relativis¬ 
tic treatment, such as will ultimately be necessary for ac¬ 
curacy. Relativistic effects will be particularly important 
for the core, in which however, compared with the crust, 
the effect of entrainment is expected to be much more 
moderate |33IIT$] at least in the outer part for which the 
basic physics is reasonably well understood, unlike the in¬ 
ner core for which many exotic possibilities have been en¬ 
visages, such that of hyperon crystallisation ESI, or of a 
LOFF phase in a color superconducting quark condensate 
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Figure 1: Intersection of rotating neutron star model with 
a Taylor cylinder, over which the superfluid angular veloc¬ 
ity must be uniform when the conditions for the Proudman 
theorem are applicable. 

2 Two fluid model for stratified 
neutron star crust 

2.1 The Lagrangian master function 

The previous treatment referred to above m used a non 
relativistic Newtonian description in terms of free neutron 
and confined baryon 4-currents rif and n " respectively 
representing a neutron superfluid, and a normal crust com¬ 
ponent consisting of ionic nuclei characterised by a con¬ 
fined baryon number A c and a positive charge number Z 
(in a degenerate electron background) that were treated 
as decoupled but interpenetrating barotropic perfect fluids 
interacting only via vortex pinning and long range gravita¬ 
tional forces, with the solid rigidity of the crust component 
taken into account by the inclusion of an extra stress force 
density. 

One of the improvements introduced here is to drop the 
barotropy restriction for the ionic constituent by allowing 
its behavior to depend on the nuclear mass number A c 
(which will tend to be larger in the deeper layers) thus 
allowing for the possibility of a stratification effect that 
may be important for stabilisation. 

Another improvement that would be needed for a more 
accurate analysis is of course the use of a general rela¬ 
tivistic treatment, of the kind that has long been com¬ 
monly used with a simple perfect fluid description, and 
that has recently been developed for use with a two con¬ 
stituent fluid description of the kind needed fore treating 
the stellar core j501f51. However for the semilocal modeli- 
sation needed for the treatment of differential rotation in 
the crust layers with which we shall be concerned here, it 


would appear that corrections from general relativity will 
be much less important (at most a few tens of per cent) 
than those due to the entrainment. A recent program of 
work EH EH EH (using a microphysical analysis combin¬ 
ing methods from nuclear and solid state physics) on the 
relative motion of unbound neutrons in the inner crust lay¬ 
ers (between the neutron drip density about 10 11 g/cm 3 
and nuclear saturation density about 10 14 g/cm 3 ) has pro¬ 
vided results EH indicating that the effect of entrainment 
is likely to be very much stronger than had previously been 
supposed. 

As in the preceding work EH, we shall therefore pro¬ 
ceed here within a non relativistic Newtonian treatment, 
but instead of using a decoupled perfect fluid description 
we shall now include allowance for the effect of entrain¬ 
ment, using a non-dissipative two constituent description 
of the generic kind already developed for use in the liquid 
core El El EH where the entrainment corrections are rel¬ 
atively moderate, though still potentially important as a 
cause for two-stream instability |Hj if the relative velocity 
becomes too large. 

For the reasons outlined in the introduction, it is both 
instructive and convenient (even though our treatment is 
non-relativistic) to use a canonical formalism of the 4- 
dimensionally covariant kind recently developed El EH 
for the purpose of facilitating the exploitation of the varia¬ 
tional formulation that is applicable to the non-dissipative 
limit. In this limit there is no resistance to relative mo¬ 
tion, so that the only direct coupling between the two 
constituents is via the entrainment effect El , whose micro¬ 
scopic origin is the Bragg scattering of dripped neutrons by 
crustal nuclei EH which is the analog in the nuclear con¬ 
text of conduction electron scattering in ordinary solids. 

It is also convenient to use a chemical constituent label x 
that runs over two values, namely x = f for the superfluid 
constituent and x = c for the crust constituent. Each 
constituent will have a baryon current density 4-vector 
having the form 

K = n x< (!) 

in which the scalars n x are the respective number densities 
of free baryons (only neutrons) and of baryons that are 
effectively confined within the atomic nuclei (namely all 
the protons, together with a roughly equal but somewhat 
larger number of confined neutrons, whose precise specifi¬ 
cation needs further clarification as discussed below) while 
the vectors u” are the corresponding 4-velocities, whose 
components are labelled by a Greek index v running over 
the values 0,1, 2, 3. The conservation of the total baryon 
number current 

< = <+<, ( 2 ) 

is concisely expressible in this 4-dimensional notation as 

= 0. (3) 
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In Aristotelian coordinates (representing the usual kind 
of 3+1 space time decomposition with respect to the rest 
frame of the star) the 4-velocity components will have the 
form = 1, = tP where the v £ are ordinary 3-velocity 

components, using a Latin index i = 1,2,3. This means 
that time component of the current density is simply equal 
to the corresponding particle number density n° = n x , 
while the space components are the those of the current 
density 3-vector n 1 . 

In addition to the six independent vectorial components 
provided by the space vectors n* and rif and the pair of 
scalars provided by the “confined” and “free” baryon num¬ 
ber densities n c and m, our model will contain a nineth 
independent variable, namely the scalar provided by the 
nuclear mass number A c , or equivalently by the number 
density m of nuclei, whose number current will be given 

by 

Tii = n i u c > ni=n c /A c . (4) 

We shall take the mass per baryon to be the same (con¬ 
sidering the electron mass and the proton neutron mass 
difference to be negligible compared with the ordinary 
proton mass) with fixed value m c = m f = m, so that 
the corresponding mass density contributions will simply 
be given by p x = mn x . The total mass density, as given 
explicitly by 

p = mn B = p { + p c , (5) 

will be important as the Poissonian source of the Newto¬ 
nian gravitational potential scalar, <f>. Within this New¬ 
tonian background field, the total force balance equation 
will take the form 

v M n = /„-pv^, (6) 

in which /„ is the non gravitational external 4-force den¬ 
sity - if any - that may be acting on the system, and T ^ is 
the stress momentum energy density tensor as constructed 
in terms of the 4-dimensional formalism recapitulated be¬ 
low. 

The dynamics of system will be governed - in the man¬ 
ner described in our preceding work mm by a La- 
grangian master function 

A = A mat (j> p , A mat = Akin + Ai n t (7) 

in which Akin is the ordinary Galilean frame dependent 
kinetic energy density, namely 

Akin = | {pfv i + p c v c 2 ) , (8) 

and Aint is some appropriately specified internal action 
density contribution that must be Galilean frame indepen¬ 
dent, which means that it can depend only on the scalar 


number densities m, n{, n c , and on the velocity difference 

K ]=vi ^ vi (9) 

so that in terms of an “entrainment density” function p ic 
and a set of chemical potentials y f , y c , x 1 its generic vari¬ 
ation will take the form 

<5 A int = Pfcbi] 8[v l ] - X l Sn{ - x c (5n c - X^m , (10) 

which gives the corresponding formula 

6A mat = nfp\Svf +n c p c i 6v^-p t Sn { -p c 6n c -p I 5ni , (11) 

in which the (Galilean frame dependent) material energy 
coefficients are given in terms of the corresponding (frame 
independent) chemical potentials by 

p { = X f - \tnvi , p c =x c -^mvc, p 1 = x* , (12) 

and the associated space momenta are given by 

p\ = mvu + —[vi \, p\=mv C i-—[vi\. (13) 

Tlf Tl c 

In the low relative velocity regimes that are relevant it 
will always be possible to take the internal energy contri¬ 
bution to have the quadratic form 

Aint = \p ic [v? - Uins , (M) 

with p[ c and the static contribution U j ns specified by ap¬ 
propriate equations of state as functions only of the scalar 
densities nj, n c , rif. It is to be remarked that to translate 
our notation to that used by Prix et al. m we would 
need to make the conversions x f i—> //, x c i—> p c and 
Aint e-i► — £ where £ is what these authors have referred 
to rather loosely as “internal energy” density, although it 
is something that should not be confused with the true 
internal energy density U i n t, which will actually be given 

b y _ . 

U int = jftcH 2 + Uins = P{ c [vf - A int . (15) 

2.2 Currents and their conjugate mo¬ 
menta in 4 dimensions 

Instead of treating quantities such as energy and momen¬ 
tum separately, as has traditionally been done in a New¬ 
tonian framework, it has been found to be technically 
advantageous, as remarked above, to unify them in a 4- 
dimensional treatment DUE! (which was originally in¬ 
spired by the corresponding relativistic theory [313). Tlie 
first step in this 4 dimensional treatment is to rewrite the 
generic variation m in terns of 4-current variations in 
the form 

dA mat = drif + p, c v d n" - p 1 dm , (16) 
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which provides a direct specification of the material 4- 
momentum covectors pr u and p' : v that play a crucial role in 
the concise canonical formalism to be described below. In 
terms of the quantities introduced above their time com¬ 
ponents will be given by 

= ~(~A + mW ), A 0 = ~{~A + M> c) > (17) 

while their space components will simply be given by 03 • 
These material 4-momenta, and p c „, give rise to cor¬ 
responding complete “free” and “confined” baryon mo¬ 
menta, 7 r f „ and 7r c „, whose time components include al¬ 
lowance for gravitational energy according to the specifi¬ 
cations 

7T o = A 0 - ™</> j tt c 0 = M c 0 - rncf ) , (18) 

while their space components are given simply by 

7T f i = At , A = A • (19) 

It is also useful to introduce the total ionic momentum 
covector tt\, which is defined so as to include allowance 
for nuclear energy by 

= A c tt c o — A , tt\ = A c w c i . (20) 

This enables us to express the complete stress energy ten¬ 
sor of the system in the very concise form 

n = n f V 1/ + <7r I I/ + $<y'; J ( 21 ) 

in which the material pressure function is given by 

^ = A mat - rif A v ~ n c A v + n i A • ( 22 ) 

It follows that the force balance equation © can be rewrit¬ 
ten in the form 

f v + f\ = f v , (23) 

in which the 4-force density / f „ acing on the “free” neu¬ 
trons is specified in terms of their vorticity tensor zUf lv by 

f v = A + A v V M n f M , = 2^7^ , (24) 

while the combined 4-force density acting on the ions and 
their “confined” baryons will be specified by the formula 

f\ = 2<V [m A v] + <7T C ^A C . (25) 

This can be decomposed in the form 

f\ = fl + niV^ 1 - S°MAK ), (26) 

in which the last term contributes only to the time com¬ 
ponent, while the only part that remains in the absence 
of stratification is the contribution of first term, which is 
given by an expression of the standard form 

/ c „ = 2<V [/1 < ] +<V M n'‘. (27) 


2.3 Equations of state and effective mass 

To complete the specification of the variables introduced 
in the previous section it is evidently necessary to pre¬ 
scribe the particular form of the material action density 
A m at, which is - for this purpose - most conveniently de- 
composible in the form 

A ma t = U dyn t^ins , (28) 

wherein U i ns is the static internal energy contribution in¬ 
troduced in GJ, which is specified by an appropriate “pri¬ 
mary” equation of state as a function of the three relevant 
number densities nt, n c ni, while the remaining dynamical 
energy contribution is given m in terms of the “normal” 
(crust) velocity 3 vector v l c and the relative current 3 vec¬ 
tor 

n l = nf[v l ], (29) 

by the formula 

U dyn = pVciVl/2 + m TliV ® + 71,71* /2K, , (30) 

in which the 1C is the “mobility coefficient”. This quantity 
1C is given by a “secondary” equation of state as a function 
of the relevant number densities ni, n c , and rif. 

In the analysis that follows a particularly important role 
will be played by the “free” neutron 3-momentum covector 
A for which, on a local “mesoscopic” scale (small com¬ 
pared with the intervortex separation) the superfluidity 
property entails the irrotationality condition = 0. 

It follows from <I3U1) that the contravariant version (as ob¬ 
tained by contracting with the ordinary Euclidean space 
metric y*- 7 ) of this (superfluid) momentum will be given 
by an expression of the form 

//* = m*t; f + (m — to*)t;* (31) 

in terms of an effective mass variable, 

771* = 71f//C , (32) 

whose difference from the ordinary baryon mass m is pro¬ 
portional to the “entrainment density” 

p {c = pi ( to * — m)/m , (33) 

which - in much of the inner crust - now seems likely m 
to be much larger than had previously been supposed. 

The reason why it is convenient for the formulation of 
the equations of state to start with the specification of 
1C rather than the effective mass 771* or the corresponding 
“entrainment density” p {c in the decomposition Hu is that 
- unlike 777* and p ic - the “mobility coefficient” 1C has the 
advantage of being physically well defined in a manner 
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that does not depend on the choice of the chemical basis 
that will be discussed in Section EZD 

Another kind of decomposition that is also independent 
of the choice of chemical basis, and that is instructive for 
the purpose of comparison with the notation used in re¬ 
lated work, including the seminal article of 0, is to intro¬ 
duce the concept of the so called “superfluid 3-velocity” 
V) s , which in stricter terminology should be referered to 
as “superfluid momentum per unit mass”, as defined by 
setting 

Mi =mVf V si = v l c + — y] . (34) 

This enables the dynamical energy density (® to be 
rewritten (without a cross term) in the seductively sim¬ 
ple form 

t/dyn = y (n s v?v sl + n N v ci v l c ) , (35) 

in which the so called “superfluid particle number den¬ 
sity” n s and the corresponding “normal particle number 
density” n N are defined by 

n s = m/C = Mf/m* n N = n B — n s . (36) 

2.4 Dynamical Equations 

Even after the primary (3 variable) equation of state func¬ 
tion t/j ns {ni, n c , rif} and the secondary (2 variable) equa¬ 
tion of state function K.{ni,n c } have been specified, the 
formulae given so far will merely provide a coherent set 
of definitions, but - except for the total baryon conser¬ 
vation law 0 - they contain no dynamical information. 
To govern the evolution of the nine independent variables 
(namely rq and the components of the 4-vectors rif and 
n") of the system we still need eight more conditions, 
which can be provided by specification of the values of 
the 4-forces }\ and /(,. 

In the absence of an external force f v in (El the sum 
of the contributions / { v and f c v would have to vanish, a 
condition that will automatically be satisfied if these con¬ 
tributions are specified just by the variational principle, 
according to which they should each vanish separately, 

fl = fl = Cl¬ 
in a dissipative model o the kind recently developed by 
iT2] (such as might be needed to allow for the mutual re¬ 
sistivity that would be present if the temperature were too 
high for superconductivity) the contributions / f „ and f c v 
could still add up to zero even though they would no longer 
vanish separately. However we wish to consider more gen¬ 
eral situations in which a non-vanishing total force density 
f v may be present for two reasons. One reason is to al¬ 
low for the “secular” (very long timescale) effect of weak 


magnetic braking of pulsar rotation. However the main 
reason for introducing the external force contribution f v 
in the present analysis is to take account of the presence, 
due to the solidity of the crust, of additional stresses that 
we wish to estimate but for which a complete treatment 
would be beyond the scope of the purely fluid model used 
here. 

In conjunction with ©, strict application of the varia¬ 
tion principle would entail as consequence that each of the 
three currents rif, n c " should be separately conserved. 
This would no doubt be a very good approximation for 
the purpose of treating high frequency oscillation modes. 
However, in the medium to long timescale processes we 
wish to consider here, whereas it will still be reasonable 
to postulate conservation of the number of ionic nuclei, 
meaning that we we shall have 

= 0, (37) 

on the other hand it may be realistic to allow for the 
possibility that - due to beta processes whereby neutrons 
are transformed into neutrons or vice versa - the “free” 
and “confined” baryon currents may not be separately 
conserved, but will evolve in such a way as to diminish 
the magnitude of the relevant chemical affinity [03, A. say, 
meaning the chemical potential difference between the two 
constituents, as measured in the relevant “normal” (non 
superfluid) rest frame, namely that of the ionic lattice with 
unit velocity 4-vector u ” , so that it will be given by 

■4 = <[/V I) W = Ml - • (38) 

As discussed in the treatment of the analogous relativis¬ 
tic model [3U (with tl\ and m 1 respectively replaced by 
entropy density and temperature) a non-dissipative model 
that is self contained - so that total force /„ vanishes - will 
be obtainable in a manner that is compatible either with 
the short term conservation conditions V !/ nf' = = 0, 

or alternatively with the long term equilibrium condition 
A = 0 that is more relevant here, provided the three force 
acting on the superfluid neutron current is specified in 
such a way that 

n/* 07 % ,nZ = 0 . (39) 

2.5 Superfluidity condition 

On a mesoscopic scale (large compared with the ionic spac¬ 
ing but small compared with the distance between quan¬ 
tised vortices) the superfluidity property of the “free” neu¬ 
trons entails that their momentum covector ir v must be 
locally proportional to the gradient of a quantum phase 
scalar, and hence that their vorticity 2-form must vanish, 
w t flv = 0 , which automatically ensures that the condition 
0| will be satisfied. 
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Figure 2: Sketch of the 2-surface swept out by a quantised 
vortex line. 

What we are concerned with here however is the macro¬ 
scopic level - meaning lengthscales large compared with 
the intervortex separation - at which there will be a non¬ 
vanishing vorticity 2-form, and a corresponding spacelike 
vorticity vector 

w x = I, (40) 

whose components in Aristotelian coordinates are simply 
w° = 0 and w 1 = e y fc V jV ik, representing the mean density 
of circulation around the quantised vortices. The superflu¬ 
idity property does however entail an algebraic restriction 
to the effect that we must have 

= 0 , (41) 

a condition need not be satisfied for a more general fluid 
motion (such as would be possible for the free neutrons in 
a very young neutron star that has not yet dropped below 
the temperature of the BCS pairing transition). This re¬ 
striction expresses the requirement that, instead of having 
matrix rank 4 as in the normal case, the antisymmetric (so 
necessarily even ranked) vorticity tensor should have 
its rank reduced to 2 in the superfluid case, a condition 
that is mathematically necessary for it to be tangential 
to a set of 2-surfaces which in this instance are swept out 
by the quantised vortex lines as illustrated on figure |5J 
This implies the existence of vorticity transport 4-vector 
field u ” subject to the usual time normalisation condition 
= 1 characterised by the defining property 

= 0 ( 42 ) 

which does not fix it completely but evidently allows a 
freedom of adjustment by addition of a spacelike vector 
field aligned with the direction of the vortex lines as given 
by w". 


The supposition that the vortices will be “pinned” in 
the sense of being transported with the ionic lattice is 
evidently expressible by adopting a dynamic equation of 
the form <o with 

< = < (43) 

which still automatically ensures that the condition in 
will be satisfied, but which will generically entail the in¬ 
volvement of a non-vanishing force density that will be 
attributable to the Magnus effect as discussed in the ap¬ 
pendix. This condition m particularly plausible as a 
description of evolution on very long timescales, for which 
the beta equilibrium condition, 

A = 0, (44) 

will also be appropriate. 

For evolution on the more moderate timescales that will 
be relevant for many astrophysical processes the situation 
is not so clear. An obvious - though not necessarily re¬ 
alistic - alternative way of satisfying ilTOli is to apply the 
strict variation principle to the effect that / should van¬ 
ish, which instead of as well as the separate “free” 
neutron conservation condition V„Tif' = 0 also provides a 
dynamical equation of the form m but with 

< = K , (45) 

and which instead of (PI) provides a separate conservation 
law 

= 0. (46) 

However the realism of such an alternative is question¬ 
able on the basis of microphysical considerations which 
suggest that both the vortex slipping implied by iTOjIi and 
the beta disequilibrium entailed by m would in prac¬ 
tise give rise to dissipative processes that would require 
a more elaborate treatment. This issue is hinted at, on 
a purely macroscopic level, by the consideration that the 
actual physical meaning of the simple set of non dissipa¬ 
tive evolution equations provided by ilTOl) and {H is open 
to question, due to its gauge dependence in the sense dis¬ 
cussed in the next section. 

2.6 Chemical basis dependence 

A problem with the force free postulate, namely the use 
of 14511 in conjunction with (Tl(l is that (unlike the use of 
<G3 in conjunction with m which will be appropriate for 
longer timescale processes) its physical meaning depends 
on how many of the neutrons are considered to be “free”. 
To be more explicit, it depends on the choice of the dimen¬ 
sionless parameter a c in the chemical base transformation 

n" = a c rip 1 ' iif = n" + (1 - a c )n p v , (47) 
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whereby the “confined” and “free” baryon currents are 
specified in terms of the more obviously well defined chem¬ 
ical basis constituted by the proton current n p " and the 
total (“confined” plus “free”) neutron current n". 

Obvious simple possibilities for the choice of such a 
“chemical gauge” include the “comprehensive” gauge, a c 
in which all the neutrons are counted as “free”, the “pair¬ 
ing gauge” a c = 2 in which the most strongly bound neu¬ 
tron states, namely those paired with proton states, are 
counted as “confined”, but all the others are considered to 
be free. A more intuitively “natural” - but more fuzzily 
defined - possibility would be to reserve the qualification 
“free” for neutrons that are actually located outside the - 
somewhat blurred - boundaries of the nuclei that contain 
the protons. However what really matters is not location 
in ordinary space but in phase space: the “operational” 
criterion for a neutron state to be effectively “free” is that 
it should have sufficient energy to overcome the barriers 
separating the ionic potential wells either classically or by 
quantum tunnelling within the relevant timescale. 

On this “operational” basis it is clear that all the neu¬ 
trons will be effectively “free” in the fluid core and that 
none of them will be effectively “free” in the outer crust 
below the “drip” density threshold, while at densities just 
above this threshold there will be a clear cut critical energy 
above which the neutrons will be able to travel over ionic 
separation distances on a microscopically short timescale 
and below which their wave functions will be exponentially 
suppressed outside the ionic wells so that they will be able 
to tunnell only on cosmologically long timescales. How¬ 
ever in the deeper and denser layers of the crust (which 
are likely to be particularly important for glitch processes) 
this “operational” distinction will no longer be so clear cut, 
because the ionic wells will get too close for the exponential 
suppression outside to be fully effective, so that there will 
be marginally bound states with intermediate penetration 
timescales that are macroscopically long but cosmologi- 
caly short. This means that the most appropriate way of 
precisely defining an “operational” basis will depend on 
the timescales involved in the astrophysical context under 
consideration. 

In view of the debatability of this question about which 
basis may be most appropriate, it is important to observe 
that the counting convention used to define the chemical 
gauge parameter a c in G3 makes no difference at all to 
the specification of the “normal” velocity v* and the rel¬ 
ative current n l defined by (El and hence that it has no 
effect on any of the terms in the original specification E1 
for the dynamical energy density. It is no less important 
to observe that uncertainties about the appropriate choice 
for a c make no difference to the specification of the “free” 
momentum covector //„, nor in consequence, to the vortic- 
ity form and the so called “superfluid 3-velocity” Vf , 


which means that the quantities appearing in the rewrit¬ 
ten formula 13511 will also be unaffected. 

The rule chosen for the specification of a c and hence of 
m will evidently affect the individual terms in the formula 
EH> for the stress energy tensor, but it can be seen that 
the total, T IJ V , will be unaffected by the value of a c so long 
as it is held fixed, and even (though this is less obvious) if 
it is allowed to vary, so long as this ratio, 

a c = n c /n p = A c ni/n p , (48) 

is postulated to depend only on the nuclear charge number 
Z, which is itself definable as the ratio 

Z = rip/rij. (49) 

It thus follows from © that provided a c (or equivalently 
A c ) is chosen as some function just of the single variable Z 7 
although the choice of this function will indeed affect the 
specification of the separate force density contributions / f „ 
and J x v . nevertheless it will not have any effect on their to¬ 
tal m ■ This means that it will not have any effect on the 
complete set of evolution equations given by specifying the 
total external force f v (e.g. by postulating that it should 
vanish) in conjunction with the ionic number conserva¬ 
tion law E3 together with the (long timescale) pinning 
and beta equilibrium conditions m and EH) . However if 
TO is replaced by TO . then the physical specification of 
the system would depend on the ansatz for the function 
a c {Zj, which for very short timescales could realistically 
be taken to be simply the “comprehensive” value a c = 1, 
corresponding just to separate conservation of protons and 
neutrons. Similarly if m is replaced by TO the result 
would again depend on the ansatz for a c {Z}, but in this 
case the appropriate choice is a question needing further 
investigation at a microscopic level. 

The strategy of the following work is not to adopt any 
particular choice of gauge but to see how far it is possible 
to draw general conclusions that will be valid whatever the 
choice of the function a c {Zj and whatever the postulate 
adopted for the specification of the vortex drift velocity 
«w 62 J>- 

3 Axisymmetric configurations 

3.1 Angular momentum convection and 
the lag formula 

As in the preceding work m on this problem, we now 
restrict our attention to configurations that are axisym¬ 
metric in the sense of being invariant with respect to the 
action of a rotation symmetry generator q v of the usual 
kind, meaning one that is spacelike, g° = 0, with space 
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Figure 3: Axisymmetry generator g l , whose magnitude is 
equal to the cylindrical radius g and the unit vector v l 
directed along the rotation axis. 


components g l given in terms of the unit vector v l along 
the relevant symetry axis, with respect to Cartesian co¬ 
ordinates r l centered on the axis, by an expression of the 
form g l = £ l]k v ] Tk as shown on figure|31 This means that 
this symmetry generator will have a scalar magnitude g 
that can be interpreted as a cylindrical radial coordinate 
(such as was denoted by vo in the preceding work by Carter 
et al. 1X9) so that it will be given in terms of the relevant 
angle of latitude 6 by g = rcos9. 

Invariance of a field with respect to such an action is 
expressible in any (Cartesian or other) coordinate sys¬ 
tem as the vanishing of its Lie derivative with respect to 
the generator. In the particular case of the superfluid 4- 
momentum covector n i v the axisymmetry condition will 
therefore be expressible as 

+ TT^V^ = 0. (50) 

In terms of the angular momentum per free (superfluid) 
particle, as defined (independently of any choice of the 
chemical gauge parameter a c ) by 

(■ = Q v 'K i v , ( 51 ) 

(so that on a microscopic scale £ will be quantised as a 
half integer multiple of the Dirac-Planck constant h) the 
axisymmetry condition ifoTitl can evidently be rewritten as 

V^=^ l/ g v . (52) 

We can use this in conjunction with m to obtain a gen¬ 
eralised Bernoulli type m conservation law to the effect 
that this superfluid angular momentum scalar £ will be 
convected by (in the sense of remaining constant along) 
the vortex flux trajectories: 

<V„£ = 0 . (53) 


We are concerned here with configurations that differ 
only very slightly from a state of circularity, meaning a 
state in which the 3-velocities are all aligned with the ax¬ 
isymmetry generator g l , so that we can write 

< = n xQ l + : Qi <j_ = 0 , (54) 

where ft x is the relevant angular velocity, and the residual 
non circular velocity contribution is supposed to be 
very small except perhaps in very brief intervals during a 
glitch. The non circular part, if any, will evidently make no 
contribution to the superfluid angular momentum scalar, 
which, according to m will be given simply by 

£ = g 2 (^m i ,fl{ + (m — TO*)ft c ^ • (55) 

Since the axisymmetry evidently implies f/V;£ = 0, the 
conservation law can be rewritten in terms of the vor¬ 
tex drift 3-velocity v ^ as a slow variation rule of the form 


( = -vJ ± V t £, (56) 

with the usual convention that a dot denotes partial dif¬ 
ferentiation with respect to time at a fixed space position. 

In the rotating neutron star configurations under con¬ 
sideration it will be justifiable to use a small perturbation 
approximation, not just for the non-circular velocity con¬ 
tributions, but also for deviations of the circular part from 
rigid rotation with a uniform fixed angular velocity value, 
f2 say. This means that for each constituent we can write 

ft x = ft + Aft x (57) 

where, like V w± the deviation Aft x is to be considered as 
relatively small. It follows that to first perturbative order 
we can rewrite the slow variation law m as 

£ = -mft< ± V^ 2 . (58) 


Using square brackets for differences between the con¬ 
stituents according to the convention 


[ft] = ft f - ft c , [ft] = ft f - ft c , (59) 

we can write l|58l) more explicitly as 
• • 

m*[ft] + mft c =- ^w_l . (60) 

Q 


It follows that, after a finite time interval during which the 
vortices have undergone a small cylindrical radial displace¬ 
ment, A g say, the ensuing (local) change in the angular 
frequency lag will be given in terms of the correspond¬ 
ing (uniform) change Aft in the observable rotation fre¬ 
quency, namely that of the solid constituent ft = ft c by 
the remarkably simple formula 


A [ft] m / Aft 

ft m* \ ft g 


(61) 


9 







In the idealised case of pinning to an absolutely rigid 
solid the displacement A g would of course vanish, so that 
A[0] would be approximately determined as a function 
just of the spherical radial coordinate r by a proportion¬ 
ality relation of the form A[f2] oc m/m*, but in practice, 
even for perfect pinning, there will be deviations from this 
as the crustal rigidity modulus is expected 021 to be rather 
low (the solid structure will be effectively “squashy”). On 
the other hand if it is supposed that pinning is entirely in¬ 
effective (since the vortices may drift by either thermal 01 
or quantum fluctuations E2) then the generalised Proud- 
man theorem derived below in Subsection roi tels us that 
the angular velocity difference [O] and hence also its 
variation A[0] will be a function not of r but of the 
cylindrical coordinate g = rcosd, so in that case (to al¬ 
low for non-uniformity of m*) the contribution from the 
displacement term in m would have to be relatively sig¬ 
nificant. 

3.2 Simplified global model 

For the purpose of drawing conclusions from available em¬ 
pirical data such as that of Lyne et al. 021, rather than 
with the locally well defined but observationally inacces¬ 
sible quantities involved in a formula such as ED it can 
be more instructive in practice to work with corresponding 
crudely defined averages, using a total angular momentum 
decomposition of the form 

J = J f + J c (62) 

in terms of “free” and “confined” contributions of the form 

JX = J n^ i g l d 3 r, (63) 

with the “free” bayon momentum given by ED and 
with the “confined” baryon momentum given by tha anal¬ 
ogous formula 

if 1 = mvl - ml [v 1 ] , , (64) 

in which the effective mass function m* for the “con¬ 
fined” baryons is determined by its “free” baryon ana¬ 
logue,analogue to* = to*, via the symmetric relation 

n c (ml — to) = rif(m* — to) . (65) 

Using an approximation in which Of is taken to be roughly 
uniform, as O c (= O) must be in any case, these contribu¬ 
tions will be given roughly by 

j f ~ i s n f + / cf o c , j c ~ / cc o c + / cf o f , (66) 

in terms of coefficients that combine to give the respective 
moments of inertia of the confined and free parts in the 
form 

f = I s + I ci , I c = I cc + I ct . (67) 


In accordance with ED the separate coefficients will be 
given, in terms of the cylindrical radius h , as the volume 
integrals 



I s 


m* m g 2 d 3 r , (68) 


while, in view of the equivalence ED , the entrainment 
coefficient will be given by 


Id 



ml) n c g 2 d 3 r 



— to*) rif g 2 d 3 r, (69) 


so the subtotals will be given by expressions of the familiar 
form 

I c = to J n c g 2 d 3 r, I f = rn j mg 2 d 3 r, (70) 

which combine to give the complete stellar moment of in¬ 
ertia as 

I = I C +I { = J pg 2 d 3 r. (71) 

In a simplified treatment it may be supposed that these 
moment of inertia coefficients remain constant (i.e. that 
the effects of a conceivable radial displacement A g in ED 
can be ignored) and that the superfluid contribution to 
the angular momentum is conserved during the interglitch 
period, A J f ~ 0, which by ED evidently implies 

I ci 

AUfCs-y^AU. (72) 

It may also be supposed that the total angular momentum 
variation is negligible, SJ = 0 during the very short time 
interval in which the star undergoes a discontinuous glitch 
transformation involving an angular momentum transfer 
given by SJ { cs —5J C so that the corresponding angular 
velocity discontinuities will be related by 

<5U f (73) 

To avoid a long term build up of too large a deviation 
of the superfluid angular velocity fif from the externally 
observable value, U c = O, the averages over many glitches 
of their total (continuous plus discontinuous) change per 
glitch should be the same, 

(Afif + Sfl { ) ~ (AO + 50) , (74) 


a condition that is equivalently expessible using the nota¬ 
tion of ED as 

(5[0])~-<A[0]}. (75) 

The immediately preceding relations can be used to 
rewrite this condition in the form 


m .. a f ) 2 

(AO) I s I ’ 


(76) 
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in which the left hand side contains only quantities that 
are directly observable, with values that have been typ¬ 
ically found ESI to be given in order of magnitude by 
(<5ft)/(Aft) ~ —10~ 2 in the case of Vela type pulsars. 
This magnitude has commonly been used as a basis for 
estimating the ratio I ! /I of the moment of inertia of the 
superfluid component involved to that of the whole star, 
on the basis of the usual supposition that there is no need 
to distinguish between the coefficients P and I s . However 
it can be seen from P|) that the relation between these co¬ 
efficients will actually be expressible in terms of a suitably 
weighted mean value m* of m* by the formula 


in such a stationary configuration will satisfy the rigidity 
condition, meaning that ft c is uniform with the externally 
observed angular velocity ft, 

Vi ft c = 0 , ft c = ft , (80) 

it can be seen that the space parts of the force density 
formulae and pot will be expressible simply as 

/ f i=n f (ViC f +^Vft f ) , (81) 

f\ = n c V i C c + n 1 V i jj l , (82) 


I s 


m* /f 

m 


(77) 


in which, according to the recent work El referred to 
above, the ratio ro*/m is likely to be substantially larger 
than unity, which leads to a correspondingly augmented 
estimate, 


P m* (5ft) 

I m (Aft) ’ 


(78) 


for the relative value of the moment of inertia P of the 
relevant superfluid component as compared with the total 
moment of inertia I of the neutron star. 


4 Stationary circularly symmetric 
configurations. 

4.1 Total stress force in the general case 


using the notation of Prix et al. m for the comoving 
particle energies, which are concisely specifiable for the 
respective “free'’ and “confined” constituents as the cor¬ 
responding diagonal components 

C { = £ { { , C c = £ c c , (83) 

of the (non-symmetric) energy matrix defined by the for¬ 
mula 

£* = -7r*< . (84) 

This is another example of a formula that becomes more 
complicated and less heuristically meaningful when in¬ 
stead of using 4 dimensional notation it is written out 
using the traditional 3 dimensional notation, in which for 
example the comoving particle energy for the “free” con¬ 
stituent will be given by the formula 

Cf = mcj) + x f - \m v? , (85) 


We now investigate the equilibrium conditions that need to 
be satisfied in the short and medium term when we neglect 
the slow secular evolution effects considered in the preced¬ 
ing section and consider the system to to be in a state that 
is exactly stationary, and with flow that is strictly circular 
in the sense that convective velocity contributions v x ' ± in 
(EJ) are absent, so that we simply have 

vi = ft f Q i , vl = ft f g\ (79) 

For a single-fluid model in such a non convective state the 
well known Proudman theorem tells us that the angular 
velocity must be uniform over each “Taylor cylinder” (as 
characterised by a fixed value of g) for the case of a perfect 
fluid that is barotropic, meaning that there is no stratifi¬ 
cation. It will be shown in the present section how this 
result can be generalised to the two-fluid context under 
consideration here. 

Using the circularity conditions C9 in conjunction with 
the axisymmetry conditions exemplified by Et , and as¬ 
suming that the solidity property of the crust will ensure 
that the “normal” flow velocity of the “confined” current 


in which, by m and El, the relevant chemical potential 
will be given by 

X f = dUi ns /dn { - ^[v] 2 dp ic /dn { , (86) 

with dp fc /dnf = ni dm^/drif + m * — m. 

Without going into such details, it can be seen that 
the gravitational contributions to these comoving particle 
energies will cancel out in their difference, which will be 
given in terms of the chemical affinity El and the angular 
velocity difference [ft] = ftf — ft by 

C c -C( = A + l\Q\. (87) 

If the crust were purely fluid, the force density /£ would 
just have to be equal in magnitude and opposite in sign 
to the force density acting on the superfluid, namely — f\. 
Detailed evaluation of the effects of solidity would require 
the use of a more elaborate model incorporating elastic 
rigidity, but as in the preceding work by Carter et al. El 
we can allow for the effect of finite rigidity by introducing 
a non vanishing stress force density /,; representing the 
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discrepancy according to the definition JED, which can be 
rewritten in the form 


n c \ p{ ) n c nt 


( 88 ) 


in order to obtain a difference on the right hand side that 
is particularly amenable to simplification using JED- In 
terms of the chemical disequilibrium force density term 
defined by 

ft =n c V l A + n 1 V l p 1 , (89) 

it can be seen that the formulae JED and JED lead to an 
expression of the remarkably simple form 

rl rf fX 

■I i _ J i _ i | 

n c nt n c 


(90) 


We thus obtain an expression of the same form 


in the direction of Vj£ which is approximately that of ViQ. 
Consequently the condition JED can be seen to provide us 
with a Proudman type theorem telling us that the energy 
parameter C f and the angular velocity difference [f2] must 
both be uniform over each of the deformed Tayor cylinders 
characterised by a fixed value of the angular momentum 
per “free” particle, l. Provided the difference [O] is small 
compared with the uniform crust angular velocity f 1 it 
can be seen from JED that these modified Taylor cylinders 
will differ only by a small deformation from the ordinary 
Taylor cylinders that are characterised by fixed values of 
the cylindrical radial coordinate g. 

It is interesting to consider what would happen if the 
force given by JEEJ) were also supposed to vanish, so that 
there would be no net force at all: 


/^ = 0^A c ViC' c + V i /i I = 0, (97) 


fi = ~fi +f x i+f h i, ( 91 ) 

Pi 

as that previously derived by Carter et al. m in a sim¬ 
plified treatment in which the effects of entrainment and 
stratification were not included, with the force density 
attributable to centrifugal buoyancy deficit given by the 
same formula as previously, namely 

f h i = n c [Q]Vi£ . (92) 

4.2 Generalised Proudman theorem for 

force-free superfluid case 

Leaving aside the pinned case for treatment in the next 
subsection, let us restrict our attention in the remainder 
of the present subsection to the case in which the neutron 
superfluid flow is postulated to be force-free, which means 
that we have 

f\ = 0. (93) 

In this case the chemical and buoyancy contributions on 

the right of JED will be the only sources for the stress force, 
which in this case will be given simply by 

f,-.n • /';• ( 94 ) 

Under these circumstances JED reduces to the form 

ViCf =0, (95) 

in which by JED we have 

£ = g 2 (jn . (96) 

Taking the exterior derivative of equation JED implies 
that gb' fe (Vj£)Vj [f2] must vanish so that Vi [f2] has to lie 


where A c is the nuclear mass number as introduced in JD. 
This would provide an analogue of the Proudman theorem, 
to the effect that the confined particle energy parameter 
C c and the nuclear potential function p l should both be 
uniform over surfaces characterised by a fixed value of A c . 

A noteworthy special case of the application of the fore¬ 
going theorems is to the scenario considered by Prix et 
al |89| in which there is no stratification, meaning that 
Vi p l = 0 and for which the differential rotation is also 
postulated to be uniform so that one has Vi [U] = 0, in 
which case it follows that in the absence of any (pinning 
or external stress) force both C c and C f would be uniform, 
as first integral constants of the motion, in the sense that 
their gradients too would vanish, Vi C c = Vi Cf = 0. 

Assuming that (as well as the gravitational background 
<f> that may have been provided by a Cowling type ap¬ 
proximation) the distributions of the angular velocities (or 
the corresponding angular momenta) and in our case also 
of the nuclear number density n\ (or the corresponding 
stratification potential p l ) have been specified in advance 
to have values resulting from the previous history of the 
system, then the specification of two more independent 
quantities such as C c and Cf will suffice to determine the 
system completely, as remarked by Prix et al. ED, who 
observed that their values as first integral constants could 
be fixed by the specification of the globally integrated to¬ 
tals of free and confined baryons. 

A force free configuration with separate conservation of 
the free and confined constituents is not however what is 
most realistic in a system that is stationary over a secu¬ 
lar timescale. The assumption of separate conservation of 
neutrons and protons is appropriate in the context of high 
frequency dynamical processes DSHD such as may give rise 
to a two stream instability (5 : . However in an effectively 
stationary context it is to be expected that transfusion 
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(by beta processes) would be able to establish the chemi¬ 
cal equilibrium condition iTHl) to the effect that the affinity 
A should vanish, which by m means that one will have 

C c = Ci + £ [fi]. (98) 

In the special case for which [fi] and hence (by the gen¬ 
eralised Proudman theorem) also C f are uniform it fol¬ 
lows, since the distribution of i will be highly non-uniform 
(roughly proportional to the squared cylindrical radius g 2 ) 
that C c will also be highly non-uniform. We thus obtain a 
scenario in which unlike that considered by Prix et al. m 
requires only a single constant of integration. This can be 
taken to be the globally integrated total baryon number, 
which suffices to determine the uniform value of Cf, from 
which C c will be obtainable unambiguously (without any 
need to specify any other arbitrary parameter) as the non- 
uniform field provided by substituting the relevant fixed 
values of [fi] and C f. A scenario of this mathematically 
well defined type differs from the kind envisaged by Prix 
et al. m in that, despite the fact that no vortex pinning 
is involved, it requires that the solid structure of the crust 
should provide a stress force density that according to E3 
will be given as a sum of stratification and buoyancy deficit 
forces by 

fi = niVi/x 1 + n c [fi] Wd . (99) 

This could only vanish if A“ 1 Vi/ 1 1 were equal in magni¬ 
tude and opposite in sign to [fi] V,£ which is mathemati¬ 
cally conceivable but physically implausible (except near 
the stellar equator) because l will depend mainly on the 
cylindrical coordinate g , whereas one would expect that 
the other quantities involved would depend mainly on the 
spherical radial coordinate r. 

4.3 Constants of integration for pinned 
superfluid case 

To deal with the case when the superfluid vortices are 
pinned to the “normal” constituent it is convenient to fo¬ 
cus attention not so much on the comoving diagonal ele¬ 
ment Cf in the energy matrix JHH) but on the cross term 

C = £ l c = Cf + [«*] ii\ , (100) 

which is interpretable as the energy per particle of the 
“free” (superfluid) neutrons with respect to the “normal” 
rest frame of the “confined” baryons. Unlike the diagonal 
elements Cf = —A v uf and C c = —n c v uf (and the other 
off-diagonal element) this particular off-diagonal element 
C = —7 fuf has the property of being invariant with re¬ 
spect to the changes of chemical basis that were discussed 
in Subsection I TGI 

The usual supposition that, instead of moving freely, as 
was supposed in the preceding subsection, the vortices will 


satisfy the pinning condition expressed by and m 
can be seen to be equivalent to the requirement that the 
“free” neutron current should be subject to a force density 
of the form 

fi = ~f\ (ioi) 

where /I is the Magnus force density reacting on the solid 
structure of the confined constituent, as given in accor¬ 
dance with the usual Joukowsky formula (see appendix) 

by 

f\=n f w\ j [v J ] , (102) 

which works out in this stationary circularly symmetric 
case as 

f 3 i = nt[ fi]V,U (103) 

It is to be noticed that, unlike the buoyancy deficit force 

density l|92|) to which it is proportional, this Magnus force 
density is invariant with respect to changes of chemical 
basis. When substituted in the general formula m for the 
total force acting on the system it provides an expression 
of the simple form 

fi = f x i~fl (104) 

By substituting the Joukowsky formula in m it can 

be seen that in terms of the crossed particle energy m 
the vortex pinning condition can be rewritten simply as 

ViC = 0, (105) 

which means that for any pinned configuration C will be a 
uniform integral constant, whose value can be fixed by the 
total globally integrated baryon number. In terms of this 
constant C the comoving energy per particle of the “free” 
constituent will be obtainable from TO> as the variable 
function of position given by 

c f = c-[n]e. (io6) 

It can be seen from ll9Sl) that the tranfusive equilibrium 
condition (TH1) will be expressible simply as 

C c = C, (107) 

which means that, when the specification of the system is 
completed by the imposition of this chemical equilibrium 
condition, then C c will also have a uniform constant value. 

The constancy of these quantities is interpretable as a 
modified Bernouilli theorem of a kind that is rather ro¬ 
bust: it is to be emphasised that whereas the constancy 
deduced by Prix et al. m for C c and C c in the unpinned 
case was dependent on very restrictive assumptions about 
the uniformity of [fi] and the negligibility of stratification, 
no such restrictions are involved in the derivation of con¬ 
stancy for Cf and C c in the pinned case. 
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In this pinned equilibrium case, the analogue of the for¬ 
mula TO for the stress force density that needs to be 
provided by the solidity of the crust will be given, as a 
difference between stratification and Magnus-Joukowsky 
contributions by the expression 

fi = niVifj 1 - n ( [n}Vi£, (108) 

which, unlike has the noteworthy property of being 
chemically invariant. 

The work in this section has used no approximations 
apart from the postulates of exact stationarity and ax- 
isymmetry, but if we make the realistic supposition that 
the configurations under consideration differ only slightly 
from a rigidly rotating reference state with fixed angular 
velocity then, to lowest order, we shall obtain the ex¬ 
pressions 

l = mClg 2 , Vi£ = 2mfl gViQ. (109) 

This dependence of the angular momentum (per superfluid 
particle) £ on the cylindrical radial coodinate g contrasts 
with the dependence mainly on the spherical radial coor¬ 
dinate r that one would expect for the number densities 
ni and rif, as well as for the angular velocity difference [O] 
(which, in view of (1611) . might be expected to be roughly 
proportional to the depth dependent ratio m/m*). This 
is why, as in the case of TO . it is physically implausible 
that the gradient of the stratification potential fi 1 could 
be adjusted so as to cancel the terms in m except near 
the stellar equator. 

5 Conclusions 

The foregoing work shows how it is possible to construct 
an extensive range of simple, explicitly integrable, axisym- 
metric stratified two-fluid neutron star models in which a 
small differential angular velocity [fl] can be chosen in ad¬ 
vance as an arbitrary (implicitly history dependent) func¬ 
tion of the spherical coordinates r and 6 in the pinned 
case, and as an arbitrary function just of the cylindrical 
coordinate g = r cos 8 in the unpinned case to which the 
generalised Proudman theorem applies. 

This idealised category of non-dissipative models in¬ 
cludes, as extreme case, the kind of force-free model con¬ 
structed by Prix et al. m on the basis of the postulate 
of separate conservation of protons and neutrons. Except 
in this special case it is shown that the equilibrium of the 
two-fluid model necessitates the presence of a balancing 
force (given by TO for unpinned configurations and by 
m for pinned configurations) that must presumably be 
provided by stress in the solid structure of the crust (whose 
detailed treatment is beyond the scope of the purely fluid 


description used here) which can build up only to a certain 
point at which the breakdown responsible for the glitch 
phenomenon will occur. 

The exceptional case in which the presence of such 
(glitch producing) stress was avoided m depended on 
the postulate of separate conservation of protons and neu¬ 
trons, a condition that would be realistic in the short 
timescale context of dynamic perturbations but not in 
the long timescale context of stationary configurations, for 
which the postulate of chemical equilibrium is advocated 
here as a more plausible idealisation. 

It is to be understood that the force density given ei¬ 
ther by (El or more likely by m will have to be bal¬ 
anced by a solid stress gradient, meaning that it will sat¬ 
isfy a relation of the form /, = V , j(S) J i in which ( S ) 3 i 
is a shear stress tensor that is the trace-free part of an 
elastic stress tensor S 3 i of the kind recently described in 
the framework of an elastic solid model that is already 
available m for the description of the outer crust, of 
which an appropriate neutron conducting generalisation 
will soon be ready E f° r application to the inner crust 
(with density above the neutron drip threshold) under con¬ 
sideration here. The idea is that this shear stress tensor 
( S)ij will be proportional to a trace-free shear strain ten¬ 
sor that can not excede a critical magnitude beyond which 
the solid structure will breakdown, and that a glitch will 
occur either when this critical value is reached or else, in 
the pinned case at what may be an earlier stage 

when the Joukowski-Magnus force density ifTMll reaches a 
vortex slippage limit at which unpinning occurs. 

For the study of processes occurring on intermediate 
timescales it will of course be necessary to include al¬ 
lowance for the dissipative processes that have been ne¬ 
glected in the idealised models considered here. As well 
as the controversial question of vortex slippage, EH ED 
ESI El EH a particularly important issue concerns the 
timescales EIE3ES1 of the weak interaction processes 
leading to chemical equilibrium. 
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6 Appendix 

6.1 Joukowsky formula for Magnus force 

The purpose of this appendix is to derive the formula used 
above for the force exerted on the fluid in the pinned case 
and to verify that it does indeed balance the collective 
effect of the Magnus forces exerted on the individual vor¬ 
tices. We shall start by evaluating the latter, using the 
generalised Joukowski formula that has recently been de¬ 
rived m for the force exerted on an isolated vortex by a 
mixture of irrotational perfect fluids. 

The present application is to the case of a vortex that 
is fixed with respect to the crust rest frame characterised 
by the flow 4 vector uf 1 of the “normal” constituent, and 
that is oriented in the direction specified by a spacelike 
unit 4-vector z+ say, for which the force per unit length 
will be given by 

T v = (110) 

where is the (chemical gauge invariant) relative flow 
current vector of the superfluid and C is its (chemical 
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gauge invariant) circulation as given by 

C f = f , (111) 

in which the integral is carried out along any closed path 
surrounding the vortex. 

Since the effect of many such vortices in closely spaced 
parallel array will be represented by a vorticity vector 

w fi = wv» ( 112 ) 

of the type introduced in CUli . whose magnitude w is the 
value of the total circulation per unit area, it follows that 
the corresponding force density (obtained as force per unit 
length per unit area) will be given by 

ft = ^T v (H3) 

so we end up with the formula 

ft = n pj £^ p<J w p u^ ( 114 ) 

in which it is to be observed that the value of the cir¬ 
culation C { round an individual vortex has cancelled out. 
(Thus it is not necessary to know that it will actually be 
given in terms of the Dirac-Planck constant by the formula 
C f = 7r H in which the usual factor 2 is missing because of 
the pairing due the fermionic nature of the neutrons.) 

To obtain the force density on the fluid we start from 
the formula ES in which, since we are considering a sta¬ 
tionary axisymmetric situation, the divergence term will 
automatically drop out leaving the expression 

ft = nfvotn ■ ( 115 ) 

which would be equivalent just to the ordinary Euler equa¬ 
tion if the force term vanished. It is evident that the pin¬ 
ning condition 

ufzutv = 0 , ( 116 ) 

can be used to rewrite (II Ml) in the form 

ft = M u f - KWu/j, = ■ ( 117 ) 

It can also be seen to follow from ifTTSli that, in terms of 
the vector introduced in (im the vorticity form will be 
expressible as 

™ tn = £^paW P U^ , ( 118 ) 

so that dm can be rewritten in the form 

ft = ~ nV ™tvL > ( 119 ) 

which evidently does indeed exactly balance E3- 


6.2 Effect of entrainment on vortex den¬ 
sity 

It is of interest to consider the effect of entrainment on 
the vorticity magnitude w which is the circulation per unit 
area round an infinitesimal circuit orthogonal to the local 
axis of rotation. Remarking that the circulation C f is sim¬ 
ply given by 2-7 t£, where i is the angular momentum per 
free particle, it can be seen from (1551) . by considering the 
circulation between nearby circuits with crylindrical ra¬ 
dius values g and g+Sg, that the local vorticity magnitude 
will be expressible in terms of the angular velocity differ¬ 
ence [f2] = fif — S2 C and the mass increment m f c = to* — m 
as 

w = = 2 mfi f + 2 m f c [Q] + . (120) 

ZTrgdg og 

Knowing that the circulation round an individual quan¬ 
tised vortex line in the fermionic superfluid condensate will 
be given by C { = ttH we see that the corresponding surface 
number density of quantised vortex lines will simply be 
given by w/C f = w/(nh). The first term in equation 11201 ) 
can be recognised as the well-known formula for a single 
superfluid. 

In practice, the extra terms proportional to [fi] and its 
gradient will be negligible for the astrophysical applica¬ 
tions under consideration, in which we shall generally have 
[fl] <C II, provided we are using an “operational” gauge in 
which magnitude of the effective mass m* will remain com¬ 
parable with the ordinary baryonic mass m. If the evalua¬ 
tion were carried out in a “comprehensive” gauge, meaning 
a chemical basis in which all the neutrons are deemed to 
be “free”, the relevant value of the effective mass would 
become extremely large near the “neutron drip” transi¬ 
tion, but the value of the (gauge dependent) difference 
[f2] would correspondingly become even smaller, since the 
product m*[f2] is chemically invariant, so that the same 
value would ultimately be obtained for the vorticity mag¬ 
nitude w itself. 
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